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Abstract 

We point out that oscillation between the two scalar tau (r) mass eigenstates can 
give rise to CP-violation if some parameters appearing in the stau or chargino/neutralino 
mass matrices are complex. If r + and f~ decay into different charginos or neutralinos, 
rate asymmetries as large as 20% are possible. If both staus decay directly into r+LSP, 
CP-violation can in principle still be observed through an energy asymmetry of the 
t decay products, but this asymmetry never exceeds the percent level. Even the rate 
asymmetries become small if the mass splitting between the two stau mass eigenstates 
is larger than 1%. 



1) Introduction 



The currently most widely studied extensions of the Standard Model (SM) involve Supersym- 
metry (SUSY). Superparticles can stabilize the gauge hierarchy and allow for the Grand 
Unification of all three gauge groups of the SM 0. Fortunately both of these statements 
remain true in (softly) broken SUSY models, which satisfy experimental constraints from the 
unsuccessful searches for superparticles, most notably at LEP and the Tevatron ||. Unfortu- 
nately neither of these arguments tells us anything about the dynamics, or even the scale, of 
supersymmetry breaking as long as (most) sparticles are not (much) heavier than 1 TeV. In 
phenomenological analyses it is therefore preferable to take as general an approach to SUSY 
breaking as possible, i.e. to simply parameterize it in terms of soft breaking operators. Many 
of these parameters can be complex, and therefore can give rise to many new CP-violating 
effects. 

In this note we study CP-violation in the scalar tau (r) sector at a future (linear) e + e~ 
collider. There are several reasons to assume that CP-violation might be much more prominent 
for f than for first or second generation sleptons. In the absence of generation mixing, all 
nontrivial phases in the lepton-slepton-chargino-neutralino part of Lagrangian of the Minimal 
Supersymmetric Standard Model (MSSM) can be shifted into terms involving lepton masses 
and Yukawa couplings. Moreover, present experimental bounds @ on the CP-odd electric 
and weak dipole moments of the r lepton are too weak to significantly constrain CP-violation 
in the f sector, in sharp contrast to the case of (s)electrons and, to a lesser extent, (s)muons. 
In addition, the short lifetime of the r lepton allows one to determine its polarization, at least 
in a statistical sense, by measuring the energy distribution of its visible decay products; this 
opens the possibility to construct new CP-odd observables. Finally, unlike r's, because of 
their small Yukawa couplings first and second generation sleptons could have masses in the 
tens of TeV range without destabilizing the gauge hierarchy Q; indeed, this is one possible 
U (although not very attractive 0) solution of the "SUSY flavor problem". 

In spite of these merits, CP-violation in f pair events can only be observable if the two 
f mass eigenstates are closely degenerate. The reason is that the two produced f sleptons, 
being scalar particles, decay independently of each other. All CP-odd variables are therefore 
proportional to the product of a CP-odd phase and a second, CP-even phase; this is analogous 
to the more familiar K° — K° and B° — B° systems. In the case at hand this second phase can 
only come from the Breit-Wigner propagators of the f sleptons|] or, more specifically, from 
the interference of two different r propagators. These interference contributions, which can 
be understood in terms of oscillation between the two f mass eigenstates, will become very 
small if the mass difference is much larger than the (average) decay width of the two f states. 

In this respect the situation is somewhat similar to the case of CP-violation through slepton 
flavor oscillation |7], [8] (see also refs.|| for earlier, related work). There are two differences, 
however. First, the two f current eigenstates do not carry different quantum numbers that 
are conserved in the SM, while sleptons of different generations obviously do; even detecting 



the presence of mixing is therefore not entirely straightforward in our case [10]. Second, the 
two t current eigenstates have different SU(2) x U(l)y quantum numbers. According to our 
current understanding of SUSY breaking there is therefore no good reason to assume that 
the soft breaking masses for these states should be close to each other. On the other hand, 
nothing forbids such an "accidental" near-degeneracy, either. It might therefore be useful to 



*In general there can also be other dispersive phases, e.g. due to loop corrections to the r decay vertices. 
However, they are too small to be of much use. 
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point out that such a degeneracy can lead to interesting new phenomena. 

The rest of this article is organized as follows. In the next Section we discuss rate asym- 
metries, which can occur if the positive and negative f decay into final states that are not 
CP-conjugates of each other. We find that under favorable circumstances these asymmetries 
could be detectable at future e + e~ colliders. In Sec. 3 we discuss the pion energy asymmetry 
that can result from f — > r — > it decays if both f decay into the same neutralino. We find 
this asymmetry to be too small to be detectable with the currently foreseen luminosity of 
next-generation colliders, at least within the MSSM. Finally, Sec. 4 is devoted to a summary 
and conclusions. 



2) Rate Asymmetries 

We begin our discussion with a treatment of rate asymmetries in processes of the kind 

e+e - r~f + - (xf)2 ( X f)t ■ (1) 

Here a, b labels the combination of a r with a neutralino x° = 1 to 4 in the MSSM) or 
a u T with a chargino xtn (m = 1 or 2). It is crucial to include all combinations t~t^ in the 
intermediate state, where i, j G {1,2} labels the f mass eigenstates. They are determined by 
the mass matrix in Tl, tr basis: 

2 _ / m 2 - L -m T (A T + //tan/3) \ 

J f \ -m T {A* + /Ttan/3) m\ R ) ' 1 > 

We have absorbed the D— term contributions to the diagonal entries, as well as the super- 
symmetric contributions +m 2 , into the definition of m 2 L and m\ R . Note that we allow the 
soft breaking parameter A T as well as the supersymmetric higgsino mass parameter fi to be 
complex. Defining the f mass eigenstates through (ri,f 2 ) = (tl,Tr) U?, with 

T j ( cos6 f e ict>f sin0 f \ 

I _^-«<^f COS^f ' ' 



we have 



4>t = — arg (A^ + /itan/3) ; (4a) 

2 2 

mi — mi 

tanfc = rf l L R , . (4b) 
m T \A T + /xtanp| 

The eigenvalues m 2 l2 are not sensitive to the phase <pf. 

Since there are two f mass eigenstates, the squared matrix element, or the cross section, 
for process ([I]) contains a total of 16 terms: 

- (e+e- - m: mt) = ^ • ^ • z*» (i, ^) {pf&a + pf<&t) 

■ A, k A,X>* f 1, ^, A" 2 f 1, 4, 4) D tk D tu (5) 



m~., m~., / \ m~., m|., 
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where e is the QED gauge coupling, s is the squared e + e~ cm. energy, and A(a, b, c) = (a + 
b — c) 2 —4ab. The first two factors under the sum describe e + e _ — > f~r + production. We have 
allowed for non-degenerate f states by introducing the average masses rrif ih = (m^ + mfj /2; 
nevertheless eq.(|5|) will not describe the (small) interference terms very well if mf 2 — is 
large. The coefficients Pf = (1 - P e -) (1 + P e +) and P|, ff = (1 + P e ~) (1 - P e +) describe the 
longitudinal polarization of the e ± beams. The effective couplings c^-, cf- are given by: 

L _ l-2siv?6 w s z 

ij tJ + 2sm 2 6 w cos 2 6 w s-M 2 ^ 1 } 

C S = 6ij ~ ' J^wfi' (6b) 

where the Zt^t^ couplings cfj are defined as 

cfj = i cos 2 f - sin 2 6W; cf 2 = ^ sin 2 6 f - sin 2 cf 2 = cfj* = sin (20 f ) . (7) 

The Ajfc factors in eq. (|^) describe the convolution of two Breit-Wigner propagators. We 
follow ref. [[F] in using the narrow width approximation, but allow for different decay widths of 
the two f mass eigenstates: 

(A lk ) = - J- -i , f 1 + ^ log ^ ; (8a) 

1 t r- r- 

*M4*) = l^-TxV' ( 8b ) 

m T ik l nk 1 i X ik 

where we have introduced the average widths T fik = (T^ + T fk ) /2 as well as the "r oscillation 
parameter" 

m f . - m fk 



x ik = \ Tk ■ (9) 



The second (logarithmic) term in eq. (pa|) has been introduced to improve agreement with a 
numerical convolution of two propagators for m? 2 > l.05m fl ; however, we will see that in this 
region of parameter space CP-violating effects are already very small. Notice that x k i = —Xi k , 
i.e. A ki = A* k . 

The last factors in eq.(|5[) describe f — » (xf) decays. In particular, Df k is given by: 

D? k = («4 - m\ a - m%) (L ta L* ka + R m R* ka ) - 2m fa m^ (L ia R* ka + R m L* ka ) . (10) 

Here, L ia and R ia are the fi(xf) a couplings for left- and right-handed f a , respectively. They 
can be computed easily from the f fx interactions listed in ref. [|TTJ, together with eqs.(|3|) and 
(£D describing f mixing; of course, care must be taken to allow for complex chargino and 
neutralino mixing matrices. Note that m/ n = Ri a = for r~ — > x v t decays. 
CP-odd rate asymmetries can be defined as 

a(xm = T-t\7t a Z (») 

cr(a b + ) + a(b a + ) 
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where we have used the short-hand notation a = (xf)a- From eq.(^) one finds 
A(XaX b ) oc AikAjt (c l3 cl t D a lk D b ; - c^D^D™ 

i,j,k,l 

= ^2 AikAji (cijC kl D^ k Djl — CjiC ik DjiDf k 

i,j,k,l 

= 2iJ2 AkAj&m (c ijC * kl D? k D%) 
i,j,k,l 

= -2 £ 3m (AkAjt) 3m (c^D^D**) , (12) 

i,j,k,l 

where we have used the abbreviation CijC kl = P^c^c^ + P^c^-c kl *. In the second term of the 
second line of eq.(|T2"D we have simply relabelled the indices as i <-> j, k ^ I, and in the third 
line we have used Cji = c*_-. In the last step we have symmetrized under i <-» k, I <-> j to show 
that the final result is real. (Recall that A ki = A* k .) 

Since there are four neutralinos and two charginos in the MSSM, one can in principle 
define 15 independent rate asymmetries. We found that for f masses within reach of a next- 
generation e + e~ collider operating at y/s ~ 500 GeV the two most readily accessible asym- 
metries, A(xiXi) an d A(xiX2)y are a l so the most promising ones. We see from eqs. (^2|) and 
( pE| ) that the rate asymmetries can only be sizable if the oscillation parameter xu is not much 
larger (nor much smaller) than 1. Since the total width of f eigenstates with mass < 200 GeV 
does not exceed 1 GeV, it is clear that a sizable effect can only be expected if the mass splitting 
is quite small. A numerical scan of the parameter space reveals that the rate asymmetries 
become maximal for rrif L = m fR , \A T \ = |/xtan/3|, and arg(A T ) ~ arg(/i) — 3.2; the last two 
conditions ensure that the two contributions to the off-diagonal entries of the f mass matrix 
cancel approximately (but not exactly). The condition for \A T \ can usually only be satisfied 
for tan/3 < 5, at least if we insist that the scalar potential should not have a deeper lying 
minimum where the f fields get non-vanishing vacuum expectation values. 

In Fig. 1 we show a scatter plot in the plane spanned by the two "effective" asymmetries 
A(XiXi) an d A(xiX2)y where we have fixed m? L = m fR = 200 GeV, picked random combina- 
tions of M 2 , \i and tan/3, and fixed the absolute value and phase of A T as described above. We 
have assumed that the SU{2) and U(l)y gaugino masses unify, which allows us to set both 
of their phases to zero without loss of generality ;f] we have checked that introducing a relative 
phase between these two soft breaking parameters does not change the result qualitatively. 
The "effective" asymmetries are defined as 



Mx a X b ) = A(XaXb) ■ (a-b+)+a(a+b-). (13) 

They directly determine the luminosity needed to experimentally detect an asymmetry with 
a significance of N a standard deviations: 

^(aM- (14) 

Un general, arg(Mi) — arg(M2), arg(/i) + arg(M2) and arg(A r ) + arg(Ma) are the rephasing invariant 
CP-violating phases relevant to our problem. In our case the phase of A T only enters through the t mixing 
phase 4>f defined in eq.(4a); however, the phase of y, also appears independently in the chargino and neutralino 
mass matrices. 
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where e is the overall efficiency. Fig. 1 shows that the anticipated luminosity of 500 GeV 
colliders, C ~ 20 to 50 fb _1 /yr, should be sufficient to probe at least some regions of parameter 
space, if the more complicated r decay modes can be reconstructed with an efficiency not far 
below the value of ~ 50% found in refs. fllOf for t + t~ XiXi final states where both r leptons 
decay hadronically. 

We observe some correlation between the two effective asymmetries if they have the same 
sign.0 The reason is that an enhancement of the x® mode in the negatively charged (t~) 
channel, or a suppression of this mode in the positively charged (r + ) channel, affects both 
asymmetries in the same way. On the other hand, the much weaker correlation between 
A(xiXi) an d A^XiX^i) f° r A(xiX%) < can be explained by the observation that both asym- 
metries would vanish in the limit of vanishing r Yukawa coupling, or vanishing gauge couplings. 
The final state chargino/neutralino therefore needs significant gaugino-higgsino mixing if the 
asymmetries are to be sizable. Since the parameters that describe the chargino mass matrix 
also appear in the neutralino mass matrix, large gaugino-higgsino mixing in one sector tends 
to lead to large mixing in the other sector as well. Finally, we note that the scenarios with the 
largest effective asymmetries tend to have a fairly small phase of the fi parameter, of order 0.1 
or less; this makes it a little easier to satisfy constraints from the electric dipole moments of 
the neutron and electron, either by choosing large values for first generation sfermion masses, 
or by tuning the phases of the relevant A parameters. 

In Fig. 2 we show the reduction of the asymmetry A(xiXi) as we deviate from near- 
degeneracy of the two t mass eigenstates. As discussed above, three conditions need to be 
satisfied for this near-degeneracy to occur, which means that one can deviate from it in three 
different directions: One can introduce some splitting between the diagonal entries of the 
mass matrix (||D, or one can vary the phase or absolute value of A T . These three directions 
are explored in the solid, short dashed and long dashed curves of Fig. 2, respectively, starting 
from a choice of parameters leading to a large asymmetry (~ 20%). Along the solid line, 
\mf R — m^ L \ ~ — m fl , while \A T \ varies by more than 100 GeV along the long dashed 
curve; the impression that A(xiXi) depends more sensitively on \A T \ than on m^ L — m^ R is 
therefore somewhat misleading. Finally, the phase of A T varies by about 7r/10 along the short 
dashed curve. Fig. 2 clearly illustrates that rate asymmetries will be too small to be detectable 
at a next-generation e + e~ collider if the two eigenvalues of the f mass matrix (|2|) differ by 
more than 1%. 



3) Energy Asymmetries 

If f + and r~ decay into charge-conjugate final states no rate asymmetry can be measured. 
This will be true in particular if tx\ is the only accessible f decay mode. In this case one 
can still construct a CP-odd observable involving the spins of the two r leptons. While these 
spins are not directly measurable, they affect the energy distributions of the visible r decay 
products. Here we study this effect using the simplest r decay mode, r ± — ■> tt ± u t , which also 
has the highest sensitivity to the r polarization. 

Since the CP-odd observable we want to study here depends on kinematical quantities, 
we need an expression for the differential cross section, rather than the total one; this also 
allows us to implement acceptance cuts that are needed to isolate f pair events from SM 

+ Since all asymmetries change sign when the sign of all CP-odd phases in the Lagrangian is flipped, we 
only show results for A(xiXi) > 0- 
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backgrounds. It can be written as: 



da(e + e~ 



T T 



(2tt) 5 2 12 s 
1 



B(t~ 



71 V Tj 



' sin 2 O f d cos f ^ A 3/2 1, —5 



mf, 



dn:-dE:-dn: + dE: + (pf 44* + Pg<%<%) ■ a^d^dj;. (is) 



eff „R „R*^ 



Here, Of is the angle of r~ with respect to the e~ beam direction in the lab frame, while 
the starred variables f2*± and E*± refer to the rest frames of f . Note that we have already 
integrated analytically over the phase space of the invisible decay products (neutrinos and 
neutralinos). This fixes the integration limits for E*± as 



E: 



tt^ I mm, max 



TTlf± 



1 + 



m 2 — m| " 

m 2 -± . 



mil V mi J \ m i± m i± 



ml m\ 



(16) 



with m f - = mr ik and m f + = m fjl in eq. (|H)|) . 

The structure of eq.fllBD is quite similar to that of eq.(||). However, the functions Df k now 
describe the entire f _ — > t~Xi ~^ ^ ^tXi decay chain, including all spin effects: 



7-J7T 



2ml 



mi — mi 



ml) RnR* kl — m T m^ (LnR* kl + RnL* kl ) 



mz 



mi — mi 



(R a R* kl - L a L* kl ) (2m fik El- - m\ k +m\- ml) 



(17) 



where Ln,Rn are the Xx^i T couplings that already appeared in eq. (|10D . Note that it is not 
sufficient to treat f — > r decays as incoherent sum of f — > tl and f — > tr decays. The 
interference between left- and right-handed r leptons gives a contribution of order m T , shown 
in the second term of eq.(|l7]); because this term involves the product of left- and right- 
handed r couplings it allows to produce a CP-odd phase purely from gauge contributions 
to the couplings L and R. In contrast, the other terms produce a CP-odd phase only if at 
least one factor of the r Yukawa coupling appears in the relevant combination of L and R 
couplings|] Of course, the expectation value of any CP-odd observable is again proportional 
to the imaginary part of -4 12. 

Unfortunately we found that the asymmetry in the tt + and tt~ energy distributions is 
always quite small. The most easily measured quantity is the total energy asymmetry, 



At 



(E w+ ) - (E n -) 
(E n+ ) + (E*-) 



where the pion energies are now taken in the laboratory frame. In our scans of parameter space 
we did not find a scenario where this asymmetry exceeds 1%. The "differential asymmetry" 



AeAE) 



da 
dE„ 4 



(E) 



da 
dE„- 



(E) 



da 
<£E_4 



(E) + Jh(E) 



(19) 



§The CP-odd contribution to eq.(|T^) is linear in the pion energies E*±. It therefore involves the product 
of the £"*-independent terms from one D— function with the term linear in E* from the other D— function. 
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can reach the 10% level in some bins, but such a large asymmetry always coincides with a 
very small differential cross section in the same bin. This is illustrated in Fig. 3, which shows 
the differential asymmetry (solid) and the differential cross section (dashed, referring to the 
scale at the right) for one of the most optimistic scenarios we found.[] In particular, we have 
assumed a 100% right-handed e~ beam here; for the given choice of parameters this increases 
the total cross section by about 30%, and increases the asymmetry by more than a factor of 
2. We have applied the f pair acceptance cuts listed in the second ref . ||T0|1 . This reduces the 



total accepted cross section by about a factor of 2, and especially depletes the region of small 
En, but has no effect on the asymmetry. Note that the differential asymmetry changes sign 
near the value of E n where the differential cross section is maximal. This should allow one 
to construct an optimized CP-odd variable, which shows better sensitivity to CP-odd phases 
than the total energy asymmetry flT8] ) does, which amounts to only 0.77% in the example 
shown. Nevertheless the significance of the observed effect cannot exceed that obtained from 
simply adding the effective asymmetry observed in each bin in quadrature. This gives a 
combined effective asymmetry Ae v in analogy to the effective rate asymmetries introduced in 
the previous section: 

A 2 En = JdEAl^E) 



+ 



dEyr+ dE. 



(20) 



Even for the optimistic case depicted in Fig. 3 Ae w only amounts to 0.017 fb 1 / 2 . This compares 
favorably to the product of the overall energy asymmetry ([18]) and the square root of the total 
cross section, which only gives 0.0032 fb 1 / 2 . However, even if we include events where only 
one of the two r leptons decays into the single pion mode, which increases the total available 
cross section by a factor of 1/B(t~ — > 7i~v T ) ~ 9 if the efficiencies for the other decay modes 
are similar to that for the 7r + 7t~ mode, we would still need to accumulate about 400 fb -1 of 
data to see an energy asymmetry at the level of one standard deviation. In contrast, the same 
choice of parameters yields an effective rate asymmetry A(x^xf) = 0.32 fb 1//2 , which would 
start to become visible after 20 fb -1 have been collected (assuming an efficiency e ~ 50%). 



4) Summary and Conclusions 

In this note we have studied CP-violating phenomena that could arise from scalar r oscillations 
at future e + e~ colliders. In Sec. 2 we found that rate asymmetries, which can occur if f has 
several different decay modes, can be sizable. However, several conditions have to be satisfied 
in order to get asymmetries that might be detectable at a next-generation collider. The most 
critical requirement is that the mass difference between the two f eigenstates should not exceed 
1%. In addition, one needs substantial higgsino/gaugino mixing in the chargino/neutralino 
sector of the theory and, of course, some significant CP-odd phases. 

Note that the first of these conditions is almost impossible to fulfill for scalar b or t quarks, 
since here the off-diagonal entries in the relevant mass matrices are (much) larger. One might 
hope that in case of b squarks at least their larger decay width could compensate for a larger 
mass splitting, if b — > g + b decays are allowed which proceed through strong interactions. 
However, in order to construct a rate asymmetry, at least one of the b squarks must undergo a 

*In order to reliably predict small asymmetries without generating hundreds of millions of events we have 
symmetrized our Monte Carlo phase space integration against exchange of ir + and ir~ momenta, including 
an exchange Of — > —Of. This ensures that the numerically computed asymmetry vanishes in the absence of 
CP-violation. 
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weak decay, with correspondingly reduced branching ratio. Moreover, determining the charges 
of the final state particles, which is crucial for constructing any CP-odd observable, is not easy 
for (s) quarks. We therefore do not expect third generation squark oscillations to be detectable 
at e + e~ colliders. 

In Sec. 3 we studied asymmetries in the pion energy distribution in events where both f 
decay as f — > tx\ — > WtXi- Unfortunately we found this asymmetry to be quite small even 
in the most optimistic case. The reason is that it vanishes in the limit where the tau mass and 
Yukawa coupling are set to zero. This is also true for the rate asymmetries discussed above; 
however, in that case substantial cancellations can and frequently do occur between different 
contributions to the relevant f- chargino/neutralino - v t /t couplings, which enhances the 
relative importance of the r Yukawa coupling, as illustrated in Fig. 1. No such enhancement 
occurs in the pion energy asymmetry. The differential energy asymmetry has to change sign 
for some value of the pion energy, since the total decay width of the r lepton is independent of 
its polarization. The structure of the decay amplitude implies that this change of sign occurs 
near the value of pion energy where the cross section is maximal; as a result, most events 
contribute only little to the total asymmetry. At least for the model with minimal particle 
content (the MSSM with general soft breaking terms) studied here, the energy asymmetry 
therefore remains too small to be detectable at next-generation e + e~ colliders. This remains 
true even if we allow for polarized beams, which can increase this asymmetry by more than a 
factor of two. On the other hand, at least in some regions of parameter space rate asymmetries 
should be detectable, and might provide us with the first evidence for CP-violation in the 
leptonic sector. 
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Figure 1: Scatter plot of the effective rate asymmetries defined in eq.(13). We have fixed 
rrif L = rrif R = 200 GeV, \A T \ = |/itan/3|, and $a = 0^ — 3.2. The phases of the gaugino masses 
have been set to zero. The remaining parameters have been varied randomly, in the ranges 
fOO GeV < |M 2 |, \n\ < 500 GeV, 1 < tan/? < 5 and -tt < 9^ < tt, where 6 A and 9^ are the 
phases of A T and fx, respectively. These results are for an e + e~ collider operating at = 500 
GeV with unpolarized beams. 
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Figure 2: Dependence of the rate asymmetry A(xiXi) defined in eq.(12) on the r 2 — f 1 mass 
difference. The point common to all three curves is defined by the parameters ^/s = 500 GeV, 
m fL = m fR = 200 GeV, \A T \ = 442.6 GeV, |M 2 | = 240.2 GeV, |/x| = 154.4 GeV, tan/5 = 2.89, 
and Oft = 0a + 3.2 = 0.19. Gaugino masses are assumed to be real and to satisfy a unification 
condition. The solid, short dashed and long dashed curves have been obtained by varying 
rrif R , 9a and \A T \, respectively, as described in the text. 
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Figure 3: The differential pion energy asymmetry defined in eq.(19) (solid), and the differential 
cross section after acceptance cuts [9] (dashed); the latter refers to the scale at the right. Note 
that both r leptons are assumed to decay into the single pion mode, which gives a combined 
branching ratio of only 1.2%; however, events where only one of the r leptons decays into 
this mode can also be used for this measurement. The values of the relevant parameters are: 
\A T \ = 454.6 GeV, |M 2 | = 270.1 GeV, |/x| = 152.3 GeV, tan/5 = 2.98, 9 A = -0.0358, and 
Q jx = 3.1642; the other parameters are as in Fig. 2, except that we have used an e~ beam with 
100% right-handed polarization. 
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